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ABSTRACT 

The VLTI instrument GRAVITY will provide very powerful astrometry by combining the light from four tele¬ 
scopes for two objects simultaneously. It will measure the angular separation between the two astronomical 
objects to a precision of 10 pas. This corresponds to a differential optical path difference (dOPD) between the 
targets of few nanometers and the paths within the interferometer have to be maintained stable to that level. 
For this purpose, the novel metrology system of GRAVITY will monitor the internal dOPDs by means of phase- 
shifting interferometry. We present the four-step phase-shifting concept of the metrology with emphasis on the 
method used for calibrating the phase shifts. The latter is based on a phase-step insensitive algorithm which 
unambiguously extracts phases in contrast to other methods that are strongly limited by non-linearities of the 
phase-shifting device. The main constraint of this algorithm is to introduce a robust ellipse fitting routine. Via 
this approach we are able to measure phase shifts in the laboratory with a typical accuracy of A/2000 or 1 nm 
of the metrology wavelength. 
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1. INTRODUCTION 

GRAVITY is a four-telescope beam combiner instrument for the Very Large Telescope Interferometer (VLTI) 
at the Paranal Observatory in Chile. The K-band light from four telescopes will be combined by GRAVITY 
for two objects simultaneously. By that the beam combiner will not only be capable of high resolution imaging 
with a beam width of 4 mas but will also provide narrow-angle astrometry to a precision of 10 pas - far beyond 
the current limits. One of the primary science goals targeted with these capabilities will be to probe physics in 
the Galactic Center close to the event horizon of the supermassive black hole. There, effects predicted by the 
theory of general relativity are expected to take place that cannot be observed in other, less extreme environments. 
Dynamical measurements taken with GRAVITY’S astrometric mode are the key therefore. The idea is to observe 
two objects, one science and one reference target, simultaneously and to measure their angular separation by 
determining their so called Differential Optical Path Difference (dOPD) as illlustrated in Figure]^ 
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Figure 1. Definition of dOPD. The differential optical path difference (dOPD) between two stars is the product As - B 
of the angular separation As between the objects and the interferometer baseline B. This quantity can be derived by 
measuring the distance d between the white-light positions of the respective interference patterns of the targets when 
internal and atmospheric optical path differences are subtracted.® 


However, measuring the angular separation between the targets to a precision of 10 pas not only requires 
to measure their dOPD to a level of a few nanometers but also all paths within the interferometer from the 
telescopes down to the instrument need to be maintained stable to that level. Here is where the metrology system 
of GRAVITY comes into play which will monitor internal dOPDs by means of phase-shifting interferometry. In 
the following, the basic working principle of the GRAVITY metrology system is presented. 


2. THE WORKING PRINCIPLE OF THE GRAVITY METROLOGY SYSTEM 


2.1 Narrow-angle astrometry 

The underlying idea of the metrology system is to trace the optical paths of the science light within the interfer¬ 
ometer by a laser beam, which is illustrated in Figure The laser beam travels the paths in opposite directions 
being injected behind the beam combination facility in GRAVITY and going all the way backwards through 
GRAVITY and the VLTI infrastructures to the telescopes. The laser light at first is split into two equally bright 
beams, each of which is injected into one of the two beam combiners of the reference and science object. Before 
the injection a phase-shifting device is integrated in one of the beams such that in the end the interference 
patterns of the two beams at the telescope pupils can be sampled by phase-shifting interferometry. 

For this reason a four-step phase-shifting algorithm is implemented in GRAVITY which applies four phase 
shifts A, B, G and D accompanied by the four corresponding intensity measurements at each individual telescope. 
This allows reconstructing the phases of the interfering laser beams. These phases correspond to the internal 
dOPDs which have to be subtracted from the respective phases of the science light beams measured at the 
GRAVITY detectors in order to obtain the true dOPD on sky between the targets. 

In more detail, a short mathematical derivation, considering only two telescopes for simplicity, can show that 
in such a design the dOPD needs to be calculated from three types of quantities, the science phases a and 
measured at the detectors A and B, the metrology phases and $2 at telescopes I and 2 as well as two terms 
related to the non-common path between metrology and science light, 5i and 52 '- 


dOPD = ^-^A - ^ ($2 - $ 1 ) - ^ (^2 - <5i) 

ZTT ZTT ZTT ZTT 


( 1 ) 


where \a and \b denote the effective wavelengths of the light of star A and B and \l stands for the wavelength 
of the metrology laser. Equivalent results are obtained when calculating the general case for four telescopes with 
six possible baselines. The difference of the non-common path terms and 82 will be calibrated by swapping 
the light of two nearby stars between the beam combiners such that the zero points of the metrology phases are 
found. The actual adjustment will not be applied online but in the data reduction. Thus, the relevant quantities 
that have to be tracked by the metrology during observations are its phases at each telescope which correspond 
to internal dOPDs of the interferometer between the beam combination and the telescopes. 

As already mentioned before, the phases of the metrology light are extracted by means of a four-step phase- 
shifting algorithm that shall be shortly introduced now. 
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Figure 2. The GRAVITY metrology system demonstrated in a schematic overview for two telescopes. A laser beam is 
launched which traces the optical paths from the GRAVITY beam combiners (BC A, BC B) backwards to the telescopes 
(Tl, T2). There, fringe patterns are formed by the interfering metrology laser beams and sampled by four photodiodes 
per telescope. The fringes can be moved across the diodes by the phase-shifting device implemented in the chain at 
a kHz-rate. Via this method the phase information can be extracted from the measured intensities at the diodes that 
reconstruct the fringes.!^ 

2.2 Four-step phase-shifting algorithm 

The measurement of the internal dOPDs in the optical paths from the VLTI to the instrument GRAVITY is 
accomplished by determining the phases of fringe patterns that form in the pupil planes. The corresponding 
intensity patterns are sampled at a few points with photodiodes by applying phase shifts to one of the two 
metrology laser beams. Simultaneously, the resulting intensities of the shifted fringes are measured at the 
photodiodes, which encode the OPD information. A minimum of three phase shifts and intensity measurements 
respectively are necessary to determine the intensity pattern and by that the fringe phase. In this respect, the 
GRAVITY metrology uses four shifts = {O, §, tt, with i G {A, B, C, D} in order to determine the intensity 
distribution unambiguously denoted as ABCD algorithm. 

According to the theory of interference of two beams with intensities Ii and I 2 , the intensities measured for 
the individual phase steps should lie on the distribution I{x, y, ai) as a function of position (x,y) on the detection 
area 

I{x,y,ai) = I1+ I2 + 2A//i/2sin($(x, 2 /) -I- at) , (2) 

where ^{x,y) denotes the fringe phase. 

These four shifts and intensities can be written as four such equations with three unknown variables Ji, I 2 
and $. Solving this system of equations for the variable of interest, the fringe phase $, leads to the simple 
relation below if the phase shifts correspond to exact multiples of |: 

$ = arctan . (3) 


In case the phase shifts ai deviate from the nominal values this formula can be generalized if the deviations 
are known, but it is foreseen to calibrate the nominal shifts to extract the metrology phases.l^ 

The uncertainties in the phase extraction of the metrology enters the total astrometric error budget of 
GRAVITY, which is shown in the next section. 




















2.3 Astrometric error budget 

As raised above, the simple four-step algorithm implemented in the metrology will be based on exact multiples 
of I being applied as phase shifts. This requires a calibration of the phase-shifting device. As a consequence, 
the precision of this calibration enters the uncertainty of the phase measurement. Following from the error 
propagation of Equation the phase uncertainty squared <5$^ amounts to: 





SP 


2/2 



-h 
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( 5 ) 


assuming that the four intensity measurements have equal uncertainties SI := SIa = SIb = Sic = SId and 
h = h = Iq- 

Thus, the phase uncertainty in radian is defined by the reciprocal of ^/2 times the signal-to-noise ratio 
Iq/SI. For the shortest baselines the astrometric error budget leaves roughly a few nanometers to the phase 
measurement accuracy within 5 minutes, such that all individual uncertainties introduced by the devices involved 
in this measurement should amount to fractions of this value including the error coming from the calibration 
of the phase shifter. Table demonstrates the influence of the phase shifter calibration errors of order A/4000, 
A/2000 and A/1000 on the astrometric error budget. The goal therefore is to calibrate the phase shifter to a 
nm-level of accuracy or even lower in order to reach the specified astrometric precision. 

Before the details of developing an appropriate calibration routine are presented, the properties of the phase- 
shifting device used for the GRAVITY metrology are summarized. 


Error term 

A/4000 

A/2000 

A/1000 

Phase shifter calibration error 

0.5 nm 

I nm 

2 nm 

Total metrology phase extraction error 

0.94 nm 

1.22 nm 

1.81 nm 

Total GRAVITY OPD error 

6.83 nm 

6.87 nm 

7.00 nm 

Astrometric error 

10.27 pas 

10.34 pas 

10.53 pas 


Table 1. Influence of the phase shifter calibration error on the astrometric error of GRAVITY. The numbers correspond 
to the observation case with four VLTI Unit Telescopes (UTs) during a period of 5 min. The phase shifter calibration 
error is one of several terms, such as the performance of the metrology laser and receivers, that enters the total metrology 
phase extraction error. The latter also adds up quadratically with other errors to a total GRAVITY OPD error term 
which can be translated into an astrometric error. A phase shifter calibration error of 2 nm would dominate the total 
phase extraction error of the metrology which then would become the second largest contribution to the astrometric error. 


3. CALIBRATION OF THE PHASE SHIFTER 

3.1 Phase shifter 

Our phase-shifting device of type MPX2000-LN-0.1 displayed in Figure is a component from Photline Tech¬ 
nologies. The phase modulation of optical signals is based on a birefringent lithium niobate crystal (LiNbOa). 
The refractive index of this material can be varied in proportion to the strength of an applied electric field. 
Thus, the phase of an electro-magnetic wave passing the crystal can be influenced by setting a voltage. On the 
surface of the crystal a waveguide is implemented by titanium diffusion. This diffusion zone leads to an increase 
of the refractive index such that total internal reflection keeps the injected light within the waveguide.® The 
manufacturer’s specifications of the device are given in Table together with the corresponding requirements 
from the GRAVITY metrology design. 

The birefringent crystal has two transmission axes with different refractive indices. Since the metrology laser 
light is polarized and guided in PM fibers, the polarization can be maintained in this respect when being aligned 
to one of these axes. Due to the presence of several connectors in the metrology fiber chain small misalignments 
can occur. For this reason, we spliced fiber polarizers to the input and output of the phase shifter to correct for 
these mismatches. 










Figure 3. Photline phase shifter of type MPX2000-LN-0.1 with input fiber to the left and output fiber to the right. The 
phase modulation is realized by setting a voltage to the radiofrequency input at the bottom of the phase shifter housing. 


Parameter 

Photline specification 

Metrology requirement 

Operating wavelength 

1900 nm - 2200 nm 

1908 nm 

Electro-optic bandwidth 

150 MHz 

4 kHz 

RF input power 

-20 V - -f 20 V 

0 V - 10 V 

Optical input power 

O.I W 

0.85 W 


Table 2. Specifications of the phase shifter. The device is designed for wavelengths between 1900 nm and 2200 nm. The 
electro-optic bandwidth of the phase shifter amounts to 150 MHz and the modulation range is specified from -20 V to 
20 V. For the GRAVITY metrology roughly 0 V to 10 V are used which correspond to a full wave for a wavelength of 
2000 nm and a bandwidth of 4 kHz is foreseen to operate system. Long-term tests showed that the injected power level 
of 0.85 W is withstood by the component although it is specified for a maximum input of 0.1 W by the manufacturer 
Photline Technologies. 

The phase-shifting concept of the metrology system used to determine internal dOPDs in the GRAVITY 
interferometer requires that the accuracy of the applied phase shifts should be of order I nm or less. In practical 
terms, this means that the translation between the voltage applied to the phase shifter and the resulting phase 
shift has to be known to that level. This relation can be measured in an interferometric test setup. In order 
to determine the relation between the applied voltage and resulting phase shift we tested two methods, a linear 
scan of a full wave and a phase-step insensitive algorithm, both of which will be presented here. 

3.2 Linear scan 

In a first simple approach, the main idea was to scan the range of a full wave by setting a linear voltage ramp 
from 0 V to 10 V in an interferometric test setup, which is drawn on the left side of Figure As in the final 
metrology design, the laser light is split into two channels of which one is feeding the phase shifter. Then the 
beams are combined again and the interferometric signal is read out with a receiver using the same type of 
photodiodes as our final design. 

We average over a few hundred linear scans as in Figure at the maximum rate of 500 Hz given by the 
response of our receiver design, to extract the intensity distribution of the interference pattern. Fitting a 
function of the following shape to the acquired data allows to model the desired relation between voltage U and 
phase shift a: 

I{U) = i + h+ h{U) + 2v^Jhl2{U) sin($ -b a{U)) . (6) 

The quantities Ii and l 2 (U) correspond to the transmission functions of the two interfering beams. The 
intensity Ii denotes the transmission of the unmodulated arm. The intensity I 2 is the transmission of the phase- 
shifted arm. It turns out that the transmission of the phase shifter depends on the applied voltage, i. e. is not 
independent of the applied phase step. We measure Ii and l 2 {U) with linear scans by disconnecting one of the 
interfering fiber chains before we do the actual fringe measurement. ^ and v are parameters which we have to 
introduce to adjust the bias and contrast of I which are not perfectly reconstructed by the previously measured 
h and hiU). 
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Figure 4. Linear scans of a full-wave modulation. The sawtooth voltage (dashed blue line) is applied to the phase shifter 
while the sinusoidal function (red) corresponds to the measured fringes at the receiver. 


If the relation a{U) of the phase shifter is linear, the intensity is a simple sine function and the four required 
voltages for the ABCD phase shifts are therefore easy to derive. However, phase-shifting devices can show non- 
linearities like the Kerr effect such that the function a{U) is of more complicated shape as in our case where 
a polynomial of 4th order is needed to fit the measured intensity distribution. This and other effects observed 
during the measurements strongly limit the achievable accuracy of the phase shifter calibration error and make 
the related analysis of the measured data complicated. 

As a summary the limiting factors in our measurements were the following: 

• The relation a{U) between the voltage U applied to the phase shifter and the resulting phase shift a is 
non-linear. 

• The measurements are influenced by laser power fluctuations of a few percent at frequencies from 0 Hz to 
50 Hz. 

• The transmission of the phase shifter depends on the applied voltage with a peak-to-peak modulation of 
around 10“^ as shown in Figure [^probably due to low finesse Fabry-Perot interference of Fresnel reflections 
in the phase shifter. 

• Fringe drifts of 0.2 nm per full wave occur due to environmental instabilities like temperature fluctuations 
and vibrations, which add up in the complete measurement of several hundred periods. 

• The birefringence of the phase shifter generates intrinsic contrast modulation due to polarization misalign¬ 
ments in the component. 

These effects need to be taken into account when modeling the fitting function used for the measured intensity 
distribution. In this approach, we were able to determine phase shifts with a typical accuracy of 2 nm correcting 
for the non-linear transmission function and the fringe drifts, however under the assumption that our model 
function 1(17) and in particular a(l7) describe reality well. Thus, our main concerns were: 

• Does our model describe reality well enough? 

• Are the four ABCD values that we determine from this calibration scheme still correlated regarding that 
they are not measured instantaneously and in particular suffer from environmental instabilities and laser 
power fluctuations? 

These results caused us to investigate other calibration methods which are less affected by non-linearities. 
We were able to overcome the main drawbacks mentioned above by using a phase-step insensitive algorithm 
discussed in the next section. 






Figure 5. Measurement of the transmission modulation in the phase shifter. The corresponding setup consists of the laser 
source feeding the phase-shifting device and a receiver measuring its output intensity while voltage ramps from -10 V to 
-I-10 V were applied. The plot shows the averaged and normalized intensity distribution versus the applied voltage. 
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Figure 6. Two calibration schemes for the phase shifter. On the left side a simple interferometric setup with the phase 
shifter in one of the two interfering beams is displayed. With linear voltage ramps fringes are scanned by at least a full 
wave and their intensity distribution is measured by a detector as a function of the applied voltage U in order to determine 
the phase shifts a that correspond to the ABCD steps. On the right, another setup is shown which includes two phase 
shifters. One phase shifter is used as a delay line which is continuously shifted by many different phase shifts /3 at a slow 
rate. For every phase shift P a phase step |qa — asl is applied to the other phase shifter by setting Ua and Ub- This 
phase step can be calibrated by measuring the corresponding intensities I a and 1b BiS a function of p. 












3.3 Phase-step insensitive algorithm 

The basic idea behind the phase-step insensitive approach is to shift by the four values instantaneously that we 
believe to be close to the required ABCD phase shifts, namely {0 V, 2 V, 4 V, 6 V} since a full wave is scanned 
by approximately 8 V for our wavelength. In this manner, we directly probe the phase shifts and the phase 
shifter non-linearities as well as environmental instabilities can be neglected, such that a fixed correlation is kept 
between the different steps or voltages applied to the component. With the phase shifter in one arm of the setup 
and a delay line in the other, as shown on the right in Figure]^ we used this fact to determine our desired ABCD 
steps with a phase-step insensitive method. The delay line can be used to slowly move the phase in order to scan 
the fringes while the phase shifter is constantly repeating approximate ABCD steps for each delay line shift. 

The correlation of the instantaneous ABCD shifts remains fixed independently of the delay line phase. In this 
manner, one obtains four intensity distributions A with i G {A, B, C, D} as a function of the delay line phase /3: 

h{l3) = + bi sin{/3 + ai) , (7) 

where the parameters correspond to the intensity offsets and bi to the amplitudes of the intensity modulation. 

These functions have the form of Cartesian coordinates of an ellipse in its parametric form. Plotting couples of 
these intensity distributions, A and Ij, against each other will result in an elliptic figure with a shape depending 
on different parameters of the system such as flux and phase shift Aij = ai — aj. The data can be fitted by the 
ellipse equation in Cartesian coordinates 


Elf + FIJj + Gif A Hh + KIj+L = Q 


( 8 ) 


with the coefficients (E, F, 
formulaP 


G, H, K, L), which allows for calculating the phase step by the following 

-F 
^^iEG 



Aij = arccos 


The parameters (E, F, G, H, K, L) can be written as 

K = {aibjY + (.cLjbif' — 2aibiajbj cos(Ay) — bfbf sin^(Ay) ( 10 ) 

E = b]/K ( 11 ) 

F =—2bibj cos{Aij) !K ( 12 ) 

G = bf/K ( 13 ) 

H — 2{ajbibj cos(Ay ) — aib'f)lK ( 14 ) 

K = 2{aibibj cos(Ay) — ajbf)lK ( 15 ) 

L = l . ( 16 ) 


L is fixed to be 1 in order to avoid that the parameters shrink to zero in the ^^-minimization that we run in 
order to hnd the parameters a^, Oj, bi, bj and Aij such that equationtends to zero. It is also possible to fit 
the coefficients (E, F, G, H, K, L) by direct analytical methods as demonstrated by Fitzgibbon et al. (1996)P 
but with severe drawbacks: 


• A complex inversion process is required to compute the parameters. 

• Perfect fringes are required, since any flux, transmission or contrast variation during the fringe scan can 
bias the result by up to several degrees. 

• A moderately high success rate and stability is achieved in simulations and real data when compared to a 
classical y;^-minimization. 




In this basic model we include the laser power fluctuations during the measurements P(t), which can be 
extracted from a Fourier transform of the ABCD modulation of our measurements. The measured transmission 
variation of the phase shifter used as delay line as a function of applied voltage lies within the frequency of the 
laser power fluctuations such that it is filtered out by P(t). Therefore, our final model is the following expression 

Ii{t) = P{t)[a^ + biSm{P{t) + ai)] . (17) 

With this function we determine the phase angles A.ab , ^bc and Acd in our experiment. In more detail, 
we apply linear voltage ramps consisting of 1000 phase steps of about 4 nm to the phase shifter which is used 
as delay line. Due to the natural fringe drift this realization of a delay line is not limited by discrete voltage 
sampling. The ABCD modulation is run at frequencies close to 1 kHz. The fringe drift typically occurs at 
0.1 nm/ms such that during an ABCD cycle the drift is smaller than 0.3 nm ~ A/6000 and has no signihcant 
influence on our measurement of the ABCD phase angles. 

We took series of measurements at different ABCD frequencies in the range between 0.5 kHz and 1.5 kHz 
and different dates few months from each other without specihc control of the environment to extract the 
phase shifter calibration error of this phase-step insensitive algorithm. Typically, the measured phase shifts are 
consistent within 0.2° or 1 nm between the series which also is the typical dispersion on each individual series 
of measurements. These results validate that the phase-step insensitive algorithm is an appropriate method to 
calibrate the ABCD phase steps for the metrology on a 1 nm-level. 

In comparison to the method of linearly scanning a full wave with the phase shifter, we not only achieve a 
precision of factor two higher, but we are also able to extract phase angles unambiguously. The advantage is 
that we do not need to model the interference fringes including the various complex effects that let the intensity 
distribution strongly differ compared to the theory of simple two-beam interference. The deficiencies of the linear 
scan method become obvious when directly comparing the results of both algorithms tested in Table[^ In Figure 
the determined phase steps Aab of one such series are shown as well as one of the respective ellipses together 
with the corresponding fit and residuals. 


Results 

Aab 

A-bc 

Acd 

Applied voltage step 

{0 V,2 V} 

{2 V,4 V} 

{4 V, 6 V} 

Phase-step insensitive algorithm 

88.7° ± 0.2° 

89.2° ±0.2° 

92.8° ± 0.2° 

Linear scan 

87.0° ± 0.4° 

86.6° ± 0.4° 

87.1° ±0.4° 


Table 3. Results of the phase-step insensitive algorithm compared to the linear scan. For the phase-step insensitive 
method the average phase shift values A as, Asc and Acd were obtained from series of measurements at different ABCD 
modulation frequencies and different dates during a time span of a few months. The absolute value of the applied voltage 
steps always is 2 V, but results in different phase angles due to the non-linear behavior of the phase shifter. Obviously, the 
linear scan shows systematically different as well as less precise results and fails in revealing the non-linearities probably 
due to flaws of the respective fit model. 


4. CONCLUSIONS 

The calibration of the metrology phase shifter presented here fulfills the specihed requirements for performing 
10 pas-level astrometry. We could demonstrate experimentally that ABCD phase shifts can be calibrated with 
an accuracy of typically 1 nm in a stable setup and environment. Despite non-linearities that occur in the 
component, the metrology hber chain or in the environment we were able to elaborate a method that is not as 
sensitive to these as simple scans of fringes in order to extract the relation between the applied voltage and the 
resulting phase shift. The phase-step insensitive algorithm is able to determine phases unambiguously when a 
robust fitting routine for ellipses is included. In our case a y^-minimization serves this purpose. 
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Figure 7. Performance of the phase-step insensitive algorithm. The top plot shows the fit of an ellipse to data acquired 
at a ABCD modulation of 700 Hz. The data are the voltages measured in Volt by the receiver and correspond to the 
intensities I a and Ib for the phases A and B as a function of delay line shifts. The fit line is displayed in red and 
the black data points are almost invisible due to the small dispersion. Larger dispersions are found when the setup is 
perturbed by vibrations, but no bias is introduced since they average out. In the figure below the htting residuals of the 
X^-minimization relative to the fringe continuum are shown as a function of time. The bottom plot displays a series of 
10 measurements with a mean phase step Aas = 88.9° ± 0.2°. 







These results were achieved under normal air pressure and room temperature. In GRAVITY the phase 
shifter is going to be operated in vacuum at I0“® mbar cooled to 240 K - conditions which might lead to a 
higher calibration accuracy provided that the instrument does not introduce higher noise levels. Currently, we 
analyze the instrument’s behavior in this respect. Furthermore, fiber differential delay lines are implemented in 
GRAVITY which can be used for the calibration of phase shifts and do not show a transmission dependency as 
the phase shifter. However, due to the transmission modulation of the phase shifter we might need to adapt the 
original phase extraction formula of Equation since the underlying fringe model of Equation does not take 
this effect into account. 

When the ABCD phase shifts will be calibrated to the specified accuracy within GRAVITY, an important step 
will be taken towards measuring the optical paths in the interferometer to a nm-level precision via phase-shifting 
interferometry and thus towards unprecedented astrometric errors at the level of 10 pas. 
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